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1.0  INTRODUCTION 

The  basic  approach  of  experimental  aerodynamics  is  to  strive  for  dynamic  similarity  in 
(usually)  subscale  laboratory  apparatus  and  extend  the  measured  data  to  full-scale,  free- 
flight  conditions  by  means  of  applicable  scaling  laws.  Unfortunately,  practical  obstacles 
sometimes  prevent  the  fulfillment  of  similarity  conditions,  and  this  is  a  matter  of  continuing 
concern.  The  two  principal  fluid  dynamic  scaling  parameters  in  this  work  are  the  Reynolds 
and  Mach  numbers,  and  it  has  proved  easier  to  duplicate  flight  Mach  numbers  than 
Reynolds  numbers.  For  most  current,  manned  aircraft,  it  is  typical  procedure  to  test  an 
aerodynamic  model,  with  its  major  features  faithfully  scaled  in  proportion  to  its  full-sized 
counterpart,  in  a  wind  tunnel  at  Mach  numbers  equal  to  those  expected  in  full-scale 
operation,  but  at  Reynolds  numbers  appreciably  less  than  the  flight  values.  Thus,  one  must 
not  only  determine  what  similarity  requirements  should  be  satisfied,  but  also  how  to  scale 
the  laboratory  results  to  compensate  for  the  deficiencies  in  dynamic  simulation. 

Geometric  scaling  of  major  features  generally  does  not  include  details  such  as  structural 
joints  and  fasteners.  A  possibly  more  serious  deficiency  in  similarity  is  the  usual  absence  of  a 
propulsion  system  when  aerodynamic  models  are  tested.  However,  these  are  not  the  subjects 
of  this  discussion. 

The  specific  problem  dealt  with  in  this  report  is  the  simulation  and  scaling  of  transonic, 
viscous  flows  when  models  of  aircraft  and  missiles  or  their  components  are  tested  in  wind 
tunnels  which  do  not  enable  data  to  be  obtained  at  Reynolds  numbers  equal  to  those  that 
will  be  attained  in  full-scale  free  flight.  This  statement  should  not  be  taken  to  imply  that  the 
duplication  of  full-scale  Reynolds  numbers  is  the  only  requirement  for  adequate  simulation. 
In  addition  to  the  Mach  number,  it  is  well  known  that  other  conditions  have  to  be  satisfied. 
Wind  tunnel  experimenters  must  be  concerned  with  accuracy  in  model  dimensions,  support- 
system  interference,  model  deflections  under  load,  free-stream  flow  angularity,  turbulence, 
flow  blockage,  and  effects  that  may  arise  from  failure  to  duplicate  the  unit  Reynolds  number 
and  wall-to-ambient-temperature  ratio.  However,  it  is  the  deficiency  in  Reynolds  numbers 
and  associated  effects  that  are  under  scrutiny  here. 

The  state-of-the-art  of  viscous,  transonic-flow  simulation  and  scaling  is  reviewed  and  the 
possibilities  for  improvement  in  testing  techniques  are  considered.  Prior  studies  have 
revealed  that  this  is  a  very  complex  problem,  especially  when  three-dimensional  (3-D)  flows 
must  be  dealt  with.  Therefore,  it  seems  best  to  briefly  review  the  background  and  summarize 
the  requirements  for  satisfactory  simulation  as  the  starting  point  before  discussing  what  may 
be  done  in  seeking  improvements.  Knowing  that  the  problem  is  severe,  one  must  aim  for  the 
establishment  of  a  heirarchy  of  requirements  so  that  compromises  may  be  identified. 
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2.0  FLOW  CONSIDERATIONS 

2.1  SHOCK  WAVE-BOUNDARY  LAYER  (SBL)  INTERACTION 

One  of  the  characteristics  of  transonic  flow  is  the  interaction  between  shock  waves  and 
boundary  layers.  Simulating  and  scaling  so  that  wind  tunnel  tests  will  lead  to  correct 
predictions  of  aircraft  performance  and  loads  requires  that  SBL  interaction  be  essentially  the 
same  on  the  scale  model  and  the  full-scale  aircraft.  Failure  to  match  flight  Reynolds 
numbers  causes  boundary  layer  characteristics  to  be  different  on  the  model  tested.  Skin 
friction,  boundary  layer  separation,  and  interaction  of  the  boundary  layer  with  shock  waves 
are  influenced  by  the  laminar  or  turbulent  character  of  the  boundary  layer  and  also  by  its 
thickness  and  details  of  its  structure.  For  a  given  body  shape,  attitude,  wall-to-ambient- 
temperature  ratio,  and  Mach  number,  these  qualities  depend  on  Reynolds  number. 
(Uniformity  and  steadiness  of  free-stream  flow,  as  well  as  freedom  from  wall  interference 
and  other  such  factors,  is  assumed  in  this  discussion  but  should  not  be  overlooked  in  specific 
cases.)  In  regard  to  transonic  flows,  most  attention  has  focused  on  boundary  layer 
separation  related  to  shock  wave  interaction,  particularly  on  airfoils  or  wings.  This  is 
because  shock  location  and  extent  of  separation  are  very  important  in  determining  pressure 
distributions  on  wings.  Other  regions  of  transonic  aircraft  may  also  be  sites  of  SBL 
interaction.  Examples  are  air  intakes,  deflected  control  surfaces,  and  internal  flows  at  off- 
design  conditions  such  as  may  occur  in  exhaust  nozzles.  However,  flow  over  wings  is  the 
principal  concern,  so  that  is  the  focus  of  most  of  this  discussion. 

A  number  of  very  informative  papers  on  this  subject  have  appeared  in  AGARD 
(Advisory  Group  for  Aerospace  Research  and  Development)  conference  proceedings  which 
will  be  referred  to  later.  Much  of  the  subsequent  discussion  will  be  based  on  the  flow 
features  sketched  in  Figs.  1,  2,  3,  and  4.  The  first  figure,  taken  from  Chang  (Ref.  1)  shows 
the  sequence  of  shock  wave  development  in  transonic  flows  over  an  airfoil  at  low  angle  of 
attack.  The  following  features  are  illustrated: 

1.  At  Mw  =  0.75  —  only  a  small  region  of  locally  supersonic  flow  exists  in  this 
particular  example,  and  no  shock  wave  forms.  The  wake  separates  near  the 
trailing  edge.  (More  thorough  discussion  of  separation  follows  in  the  next 
paragraph.) 

2.  At  =  0.81  •—  a  compression  shock  is  seen  on  the  upper  surface,  and  the 
effect  is  apparent  in  the  velocity  distribution.  Downstream  of  the  shock,  the 
flow  continues  subsonically  to  the  trailing  edge,  and  separation  occurs.  It 
should  be  noted  that  lambda  shocks  may  be  formed,  usually  being  associated 
with  laminar  boundary  layers  and  separation  at  the  shocks.  Typically, 
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multiple  lambda  shocks  merge  into  a  single  lambda  shock  as  Mach  number 
increases.  When  the  Reynolds  number  is  increased  and  the  boundary  layer 
becomes  turbulent,  the  oblique  leg  of  the  lambda  shock  disappears,  leaving  a 
near-normal  shock  as  shown  in  Fig.  1. 

3.  At  M  =  0.89  —  a  more  extensive  area  of  supersonic  flow  exists,  and  the 

00 

upper-surface  shock  is  located  farther  aft.  A  lower-surface  shock  is  located 
near  the  trailing  edge. 

4.  At  M  =  0.98  and  1.4  —  on  both  upper  and  lower  surfaces  the  local  flow  is 

CO 

mainly  supersonic. 

When  the  strength  of  the  shock  increases  enough  to  cause  the  boundary  layer  to  separate, 
the  extent  of  separation  can  be  sufficient  to  influence  the  pressure  at  the  trailing  edge.  This 
affects  the  airfoil  circulation,  its  pressure  distribution,  and  its  overall  performance.  Pearcey 
et  al.  (Ref.  2),  have  described  this  situation  in  detail.  Figures  2  and  3  are  reproduced  from 
their  paper  to  make  later  references  to  their  work  easier. 

Pearcey  et  al.  distinguished  between  two  flow  models,  A  and  B.  The  former  identifies  an 
SBL  interaction  wherein  the  only  separation  occurs  at  the  foot  of  the  shock  wave  and  the 
separation  bubble  grows  from  the  shock  toward  the  trailing  edge  as  shock  strength  increases 
(see  Fig.  2a).  Model  B  flow  is  that  wherein  separation  is  incipient  or  actually  present  at  the 
trailing  edge  (see  Fig.  2b).  If  the  separation  bubble  of  a  Model  A  flow  grows  rearward  all  of 
the  way  to  the  trailing  edge,  it  becomes  a  form  of  Model  B  flow.  Model  B  includes  other 
combinations  of  shocks  and  separation,  but  it  always  is  characterized  by  trailing-edge 
separation  (see  Fig.  3).  The  boundary  layer  is  assumed  to  be  turbulent,  and  flow  over  airfoils 
and  wings  is  considered  in  the  discussion  of  these  flow  models  in  Ref.  2. 

According  to  Pearcey  et  al..  Model  A  flows  are  better  simulated  in  wind  tunnel  testing 
because  bubble  growth  spreading  from  the  foot  of  the  shock  is  not  particularly  sensitive  to 
thickness  or  profile  of  the  boundary  layer,  provided  that  it  is  turbulent  at  separation. 

Osborne  and  Pearcey  (Ref.  3)  later  described  another  variant  of  flow  Model  B.  This  case 
is  one  in  which  a  small  region  of  supersonic  flow  near  the  leading  edge  induces  a  shock  which 
appreciably  modifies  development  of  rear  separation.  Figure  4,  taken  from  Ref.  3,  illustrates 
this  phenomenon  which  may  arise  when  free-stream  Mach  number  is  relatively  low,  but 
angle  of  attack  is  high.  The  influence  of  transition  fixing  when  dealing  with  this  type  of  flow 
was  discussed  in  the  referenced  paper,  but  no  easy  testing  procedure  was  recommended. 

It  should  be  noted  that  most  published  data  and  discussions  pertain  to  the  type  of  airfoils 
for  which  upper-surface  flow  receives  most  attention.  The  focus  on  the  upper  surface  is  also 
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apparent  in  this  report,  but  it  is  obvious  that  the  same  considerations  may  apply  to  lower 
surfaces  when  shocks  and  separations  become  important  there.  For  example,  supercritical 
airfoil  lower  surfaces  in  some  circumstances  deserve  the  level  of  attention  previously  devoted 
only  to  upper  surfaces  of  older  airfoil  designs.  An  informative  discussion  of  flow  over 
supercritical  airfoils  and  wings  has  been  presented  by  Yoshihara  and  Zonars  (Ref.  4). 


Figure  1.  Transonic  flow  over  an  airfoil  at  various  Mach  numbers  for  a 
2-deg  angle  of  attack.  From  Chang  (Ref.  1) 
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b.  Model  B  (i)  at  early  stage  with  localised 

shock  bubble  and  rear  separation 
(11)  at  later  stage  with  one  separa¬ 
tion  extending  from  shock  to  trail¬ 
ing  edge 


Figure  2.  Some  details  (schematic)  of  the  region  of  interaction 
between  turbulent  boundary  layer  and  shock  wave. 
From  Pearcey  et  al.  (Ref.  2) 


Green  (Ref.  5)  surveyed  the  work  on  SBL  interactions  published  prior  to  1969  and 
concluded  that  the  influence  of  Reynolds  number  (Re)  on  pressure  rise  caused  by  a  “free 
interaction”  was  unclear,  particularly  for  very  high  Re.  He  stated  that  the  then  common 
assumption  of  CpRe,1/10  =  constant  at  given  Mach  number  could  lead  to  serious  errors  if 
applied  over  a  wide  range  of  Re*.  Green  further  suggested  the  interim  adoption  of  the 
Reshotko  and  Tucker  (Ref.  6)  estimate  of  shock  pressure  rise  based  upon  the  relation  (for 
turbulent  flow), 

M2/Mi  =  0.762  (1) 


where  Mi,  M2  =  local  Mach  numbers  upstream  and  downstream  of  the  shock,  respectively. 
The  pressure  ratio  following  from  Eq.  (1)  is  determined  by  reference  to  oblique  shock  tables 
for  the  initial  Mach  number,  Mi.  Although  this  procedure  implies  negligible  effect  of 
Reynolds  number,  it  should  be  noted  that  turbulent  boundary  layer  properties  enter  into  the 
derivation  of  Eq.  (1),  and  it  is  only  by  assuming  constant  values  of  the  transformed  shape 
parameter,  H,  before  and  after  the  pressure  jump  that  dependence  on  Rex  is  avoided.  Green 
suggested  that  H  should  vary  with  Reynolds  number  but  that  the  functional  relationship 
was  not  clear  at  the  time  of  his  1970  publication. 
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Figure  3.  Flow  models  for  the  initial  development  of  shock-induced 
separation  (schematic);  turbulent  boundary  layers 
(Successive  stages  In  a  Mach  number  or  incidence 
increase)-  From  Pearcey  et  al.  (Ref.  2) 
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Figure  4.  Flow  Model  B  in  (he  particular  form  under  consideration. 
From  Osborne  and  Pearcey  (Ref.  3) 


When  Green  returned  to  this  subject  in  his  1971  paper  (Ref.  7),  he  presented  the 
qualitative  picture  in  Fig.  5,  which  is  based  on  fragmentary  data  from  several  different 
investigators.  Then,  further  to  the  question  of  shock-induced  separations.  Green  made  the 
point  that  values  of  Re0  at  a  shock  at,  say,  40-percent  chord  on  a  model  and  a  full-scale  wing 
often  can  lie  on  either  side  of  the  transition  between  “low”  and  “high”  Refl.  Green  made  it 
clear  that  Fig.  5  was  somewhat  conjectural,  and  systematic  data  were  lacking  at  the  time  of 
writing.  Since  that  date,  Settles  et  al.  (Ref.  8)  have  presented  data  that  indicate  no  Reynolds 
number  effect  on  two-dimensional  (2-D)  ramp  angles  for  incipient  separation  when  Rej  > 
105  (at  =  2.9).  The  more  recent  data  do,  however,  show  an  influence  of  Rej  on  extent  of 
upstream  influence  represented  as  Ax/Sp  For  constant  intensity  and  ramp  angle,  raising 
Rej  lowered  Ax/fy.  Assuming  Ax/5]  to  represent  scale  of  separation,  this  appears  to  confirm 
the  trend  of  curve  (b)  at  larger  Re$  in  Fig.  5. 

The  surge  in  computational  fluid  dynamics  (CFD)  during  the  past  decade  has  produced  a 
number  of  papers  on  viscous-inviscid  interactions.  Good  examples  include  the  review  by 
Melnik  (Ref.  9)  and  other  work  reported  in  the  same  conference  proceedings.  Melnik  has 
presented  CFD  solutions  for  “weak”  shocks  that  depend  on  Cole’s  wake  parameter  and  a 
normalized  shock  strength  parameter, 

Xt  =  (Mi  -  l)/e  (2) 
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where  e  =  u./Ue  and  u*  = 
parameter, 

K, 


friction  velocity  =  0  (in  Re)  - 1,  or  a  viscous,  transonic  similarity 


=  Xt/[(7  +  DMi]  = 


Mj-  1 
(7  +  l)M^e 


(3) 


where  t  =  (cn/2)1/2.  The  weak  shock  limit  is  defined  as  =  0(1),  but  results  are  given  for 
higher  xi  values.  These  results  display  a  Reynolds  number  or  scale  influence.  Low  Xt» 
corresponding  to  larger  Cfi  at  fixed  Mx  (lower  Re)  and  all-turbulent  flow,  leads  to  greater 
calculated  upstream  influence. 


Figure  5.  Influence  of  Reynolds  number  on  the  properties 
of  shock-induced  separation  at  Mach  1.5 
(conjectural).  From  Green  (Ref.  7) 


2.2  THE  CASE  OF  SWEPT  WINGS 

The  problem  already  seems  severe  enough  without  considering  the  added  complications 
presented  by  swept,  finite  wings.  In  addition  to  Refs.  5  and  7,  the  SBL  interaction  in  3-D 
wing  flows  has  been  reviewed  by  Hall  (Ref.  10)  and  Haines  (Ref.  1 1),  to  name  two  notable 
examples. 
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The  formidable  challenge  of  3-D  transonic  flows  is  apparent  when  the  basic  shock  wave 
pattern  of  Fig.  6  is  seen  to  lead  to  the  complexity  of  Figs.  7  and  8  which  are  taken  from  Ref. 
5,  and  Fig.  9,  from  Ref.  10.  Leading-edge  separation  and  vortex  formation  (Fig.  10)  is 
equally  daunting.  These  sketches  and  the  spatial  variation  that  may  exist  on  a  wing  should  be 
kept  in  mind  when  the  application  of  boundary  layer  trips  or  other  controls  is  discussed. 


Figure  6.  Basic  shock  pattern.  From  Hall  (Ref.  10) 

2.3  THE  UNIT  REYNOLDS  NUMBER 

It  is  important  to  note  that  a  form  of  unit  Reynolds  number  (U/v)  effect  has  been  found 
in  SBL-interaction  experiments  wherein  that  quantity  and  wetted  length  or  boundary  layer 
thickness  were  varied  independently.  A  constant  Reynolds  number  based  on  a  wetted  length 
still  leaves  boundary  layer  thickness  as  a  variable  influencing  separation  pressure  rise  and 
extent  of  separation.  Stanewsky  and  Little  (Ref.  12)  presented  data  on  2-D  airfoils  showing 
effects  of  unit  Reynolds  number  and  mode  of  transition  (natural  and  tripped)  on  pressure 
rise,  separation  bubble,  and  extent  of  rear  separation.  They  reported  that,  for 
incompressible  turbulent  boundary  layers, 

If  unit  Reynolds  number  is  held  constant  and  5  is  increased,  Cps  (separation 
pressure  coefficient)  will  decrease  with  increasing  Re*  (Re  based  on  boundary 
layer  thickness  upstream  of  shock).  On  the  other  hand,  if  one  holds  plate  length, 
for  instance,  constant  and  increases  unit  Reynolds  number,  Cp$  will  increase. 
(Halligan,  internal  Lockheed-Georgia  Company  report  given  as  source.) 
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However,  the  variation  of  Cps  with  Rea  shown  in  the  analytical  results  is  modest,  being 
on  the  order  of  10  percent  for  each  factor  of  10  in  Reg.  Reference  12  also  includes  data 
showing  that  extent  of  rear  separation  in  compressible  flow  over  an  airfoil  decreased  with 
increasing  unit  Reynolds  number  and  increased  with  increasing  boundary  layer  thickness. 


a.  Attached  flow 


b.  Separated  flow 

Figure  7.  Surface  flow  beneath  swept  shocks.  From  Green  (Ref.  5) 


Settles  et  al.  (Ref.  1 3)  have  found  that  2-D,  turbulent,  compressible  SBL  interaction  at  a 
compression  comer  has  an  upstream  influence  distance,  Lm,  that  scales  according  to  the 
relation. 
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Lm(U/ v)il/3  5,-2/3  =  (Lta/a1)(U«,/»)1>/3  =  const  (4) 

where  is  the  thickness  of  the  undisturbed  boundary  layer  at  the  comer  for  fixed  Mi  and 
comer-ramp  angle.  Based  on  data  and  dimensional  analysis,  the  same  investigators  propose 
another,  more  complex  relation  for  3-D  flows, 

Lm(U/*0i1/3  fi,'2'3  =  f[z(U/*-),"3  Sr2'3.  or,X,Mi]  (5) 

where  5|  is  the  same  as  $1  except  it  now  may  have  a  spanwise  variation.  The  spanwise 
coordinate  is  z;  X  is  the  sweep  angle  of  the  ramp;  and  a  is  the  ramp  angle. 


a.  Shock  parallel  to  leading  edge  on  infinite 
yawed  surface 


c.  Tapering  swept  wing 

Figure  8.  Separation  beneath  swept  shocks.  From  Green  (Ref.  5) 
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In  the  simpler  case  of  2-D  flow,  the  Ref.  13  results  show  that  upstream  influence,  i.e., 
beginning  of  separation,  scales  as  Si273  if  (U/r)i  is  held  constant,  but  Lm  will  decrease  as 
(U/v)i,/3  increases  while  5)  or  wetted  length  remains  fixed. 


Figure  9.  Flow  pattern  for  swept,  short  bubble  and  trailing-edge 
separation.  From  Hall  (Ref.  10) 


Figure  10.  Flow  pattern  for  vortex  separation.  From  Hall  (Ref.  10) 


2.4  LOCATING  THE  TERMINAL  SHOCK  WAVE 

Duplication  of  the  full-scale  shock  position  relative  to  the  chord  or  other  relevant  length 
is  a  primary  criterion  for  successful  simulation.  It  is  not  the  only  criterion,  but  it  clearly  is 
one  of  the  foremost  requirements.  If  one  is  to  do  that,  it  is  necessary  to  have  some  method 
for  predicting  where  the  shock  will  be  under  flight  conditions.  Therefore,  it  is  fortunate  that 
substantial  progress  in  that  area  has  been  made.  Cahill  and  Connor  (Ref.  14)  and,  more 
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recently,  Khan  and  Cahill  {Ref.  15)  have  presented  a  method  to  predict  the  effects  of 
changes  in  Reynolds  number  on  pressure  distributions,  once  the  distributions  are  measured 
at  lower  Reynolds  numbers. 


In  Refs.  14  and  15,  trailing-edge  pressure  coefficient  (Cpte)  is  used  as  an  indicator  of 
significant  trailing-edge  separation.  Such  separation  typically  causes  the  shock  to  move 
forward.  For  a  given  Mach  number,  a  single  curve  of  Cpte  versus  shock  location,  in  percent 
of  chord,  (Xcsh)  is  shown  to  correlate  a  large  amount  of  tunnel  and  flight  data.  These 
correlations  are,  of  course,  very  sensitive  to  Mach  number  when  that  factor  is  near  unity. 
The  parameter  used  in  Ref.  15  is  given  in  Eq.  (6).  [See  also  Ref.  9  and  Eq.  (3)  of  this  report.! 


(Me2  -  .1) 

(-Y  +  1)  e  Me2 


(6) 


where  e  =  (cf/2)l/2,  and  Me  is  the  local  Mach  number  normal  to  a  constant  percent-of-chord 
line.  Both  Me  and  Cf  are  the  values  immediately  upstream  of  the  terminal  shock. 


It  is  necessary  that  data  be  available  for  use  in  drawing  the  xCSh  versus  Cpte  curve  over 
some  range  of  Re  for  a  given  airfoil  or  wing  station  and  Mach  number .  The  parameter,  K,  is 
calculated  from  the  mearsured  airfoil  pressure  distributions  at  various  Re,  supplemented  by 
an  appropriate  skin- friction  computation  procedure.  These  results  enable  plotting  a  single 
curve  of  Cp,e  versus  K,  which  will  resemble  the  example  shown  in  Fig.  1 1.  The  break-point 
coordinates  (K0,  Cpteo)  reflect  the  Reynolds  number  influence.  Thus,  Ref.  15  contains 
curves  of  AK0  and  ACpteo  versus  Rec,  where  the  latter  is  based  on  local  wing  chord  and  the 
A*s  indicate  change  from  the  values  when  Rec  =  107.  The  experimentally  and 
computationally  determined  curve  of  CpLe  versus  K  is  shifted  to  correspond  to  higher  Rec  by 
replotting  in  (Cple  -  Cpteo),  (K  -  K0)  coordinates  and  using  the  AK<j  and  ACpteo 
correlations  of  Ref.  15  to  predict  the  Reynolds- number  effect.  The  Cp  at  the  beginning  of 
the  shock  is  calculated  using  K,  M^,  and  Cf,  with  the  stipulation  that  the  local  Mach  number 
normal  to  the  element  lines  of  the  wing  planform  at  the  foot  of  the  shock  is  unity. 


The  procedure  is  inapplicable  in  cases  where  a  local  separation  bubble  occurs  at  the 
shock,  or  low-speed  stall  occurs.  The  latter  situation  is  identifiable  with  a  smoothed-out 
compression  originating  near  the  leading  edge  and  no  evidence  of  a  discrete  shock.  The 
correlation  is  restricted  to  moderate  angles  of  attack.  In  those  cases  where  separation 
bubbles  lie  at  the  foot  of  the  shock,  significant  scale  effects  are  not  anticipated  if  a  turbulent 
boundary  layer  exists  at  the  shock.  Although  the  referenced  work  does  not  cover  every 
condition  encountered  in  simulating  full-scale  flows,  it  does  provide  an  approach  for  a  large 
class  of  flows.  Time  will  tell  if  other  users  find  the  method  to  hold  up  when  still  more 
applications  are  tried. 
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Figure  11.  Trailing-edge  pressure  coefficient  as  a  function 
of  correlation  parameter  K.  From  Khan  and 
Cahill  (Ref.  15) 

2.5  VORTEX  SEPARATION 

It  has  been  shown  that  a  swept  wing  may  experience  separation  in  the  form  of  swept, 
short  bubbles,  trailing-edge  separation,  and  vortex  separation.  Combinations  or  multiple 
examples  of  these  phenomena  may  be  present  at  one  time.  All  of  the  discussion  thus  far  has 
concerned  the  first  two  forms  because  less  has  been  written  about  scale  effect  on  vortex 
separation.  As  Hall  has  remarked  (Ref.  10),  there  will  be  some  influence  of  Reynolds 
number  whenever  there  is  flow  separation,  and  the  amplifying  three-dimensionality  as 
separation  spreads  on  swept  wings  poses  an  especially  tricky  simulation  problem. 

Leading-edge  separation  and  vortex  formation  may  be  expected  for  thin,  swept  wings  at 
rather  low  angles  of  attack.  Depending  on  the  local  Mach  number,  there  may  also  be  shock 
waves  in  that  region.  Thus,  complicated  flows  and  potentially  significant  effects  involving 
vortices,  vortex-shock  interactions,  and  other  forms  of  flow  separation  may  be  encountered. 
An  appreciation  of  this  type  of  flow  may  be  gained  by  examining  Figs.  6,  7,  and  8  in  the 
paper  by  Moss  (Ref.  16). 

When  tests  are  conducted  in  wind  tunnels  at  relatively  low  Reynolds  numbers,  the 
boundary  layer  may  be  laminar  where  the  vortex  forms  near  the  leading  edge.  If  the  full- 
scale  boundary  layer  is  turbulent  in  the  same  region,  differences  between  the  wind  tunnel 
and  full-scale  pressure  distributions  can  be  expected.  The  overall  effect  of  this  will  vary  with 
the  particular  circumstances,  and  it  cannot  be  assumed  that  the  application  of  boundary 
layer  trips  offers  a  satisfactory  solution.  The  thinness  of  the  boundary  layer  near  the  leading 
edge  makes  it  likely  that  the  use  of  conventional  2-  or  3-D  roughness  for  tripping  will  cause 
too  much  thickening  of  the  boundary  layer  and  interact  too  strongly  with  the  vortex  and 
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shock  when  those  are  present.  A  trip  suitable  for  attachment  to  leading  edges  without 
significant  thickening  of  the  boundary  layer  obviously  would  be  valuable.  Some  suggest  that 
trips  should  be  at  least  15  percent  of  chord  upstream  of  a  shock  wave. 

Rogers  and  Hall  (Ref.  17)  and  Squire  (Ref.  18),  for  example,  show  that  separation  and 
vortex  formation  on  wings  with  highly  swept  leading  edges  occur  at  very  low  angles  of 
attack,  particularly  near  the  tip.  Figure  12  from  Ref.  18  illustrates  this  in  the  case  of  a  delta 
wing  with  aspect  ratio  of  4/3  and  Reynolds  number  based  on  root  chord  of  107.  Kuchemann 
(Ref.  19)  has  discussed  vortex  and  separation  formation  generally.  He  has  stated  that 
Reynolds  number  cannot  affect  the  general  character  of  the  various  vortex  patterns  on 
swept-back  wings  but  that  it  will  greatly  influence  the  changeover  from  one  type  of  flow  to 
another,  as  well  as  shape  and  stability  of  the  vortex  sheets.  In  a  valuable  review  of 
separation.  Smith  (Ref.  20)  has  written  that, 

If  the  separation  takes  place  from  a  salient  edge,  its  occurrence  is  independent  of 
the  boundary  layer  structure,  at  least  for  high  Reynolds  number.  A  shear  layer 
forms  at  the  salient  edge,  and  again  for  high  Reynolds  number,  this  can  be 
modeled  by  a  vortex  sheet.  The  role  of  viscosity  is  to  produce  the  separation;  its 
location  is  determined  geometrically.  Inviscid  models  can  then  be  constructed. 

Hummel  (Ref.  21)  has  presented  data  for  delta  wings  with  sharp  leading  edges  in  subsonic 
flow  showing  that  the  pressure  distributions  were  different  for  laminar  and  turbulent 
boundary  layers.  However,  the  integrated  pressure  distributions  were  about  equal,  and  he 
concluded  that  the  overall  characteristics  of  slender  delta  wings  are  independent  of  Reynolds 
number. 

Polhamus  and  Gloss  (Ref.  22)  stated  that. 

The  origin  of  the  primary  vortex  feeding  sheet  and  its  feeding  rate  are  sensitive 
to  Reynolds  number,  especially  for  the  rounded  leading-edge  case  and  for  the 
secondary  vortex  regardless  of  the  leading-edge  condition.  Both  the  symmetric 
and  asymmetric  vortex  breakdown  characteristics,  which  play  an  important  role 
in  high  angle-of-attack  stability,  are  Reynolds  number  dependent. 

They  also  remarked  that  the  vortex  feeding  sheet  remains  fixed  at  sharp  leading  edges,  and 
therefore,  the  major  features  of  the  flow  are  not  critically  dependent  on  Reynolds  number  in 
that  case.  They  presented  the  data  in  Fig.  13  to  illustrate  Reynolds  number  effect  on  a  wing. 
Figure  14  is  reproduced  from  Ref.  23  to  clarify  the  primary  and  secondary  vortices  referred 
to  by  Polhamus  and  Gloss. 

In  Fig.  13  it  will  be  noted  that,  although  the  Reynolds  number  effect  is  concentrated  over 
the  forward  16  percent  of  the  wing,  the  change  in  pressure  distribution  is  marked  when  the 
vortex  is  affected. 
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Figure  13.  Effect  of  Reynolds  number  on  vortex-flow  development  —  F-lll 
TACT  flight  experiment;  a  -  6  deg,  MB  =  0.6. 

From  Polhamus  and  Gloss  (Ref.  22) 


/ 
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Figure  14.  Classical  leading-edge  vortex  characteristics. 
From  Miller  and  Wood  (Ref.  23) 
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Miller  and  Wood  (Ref.  23)  present  data  on  surface  flow  patterns  and  pressure 
distributions  near  the  trailing  edges  of  delta  planforms  in  supersonic  flows  for  various  vortex 
conditions.  The  interaction  between  vortices  and  shock  waves  in  supersonic  flow  is 
illustrated,  and  the  areas  wherein  the  different  forms  of  vortex-shock  interaction  which  may 
be  encountered  are  mapped  in  terms  of  angle  of  attack  and  Mach  number  normal  to  the 
leading  edge.  Reynolds  number  variations  apparently  were  not  investigated. 

Figure  15,  after  Waggoner  (Ref.  24),  illustrates  again  that  an  influence  of  Reynolds 
number  can  be  significant  when  vortical  flows  are  involved.  It  also  shows  that  very  large 
ranges  of  Reynolds  number  may  be  required  for  investigation  of  this  type  of  flow. 

The  seeming  contradictions  in  these  reports  can  only  be  partially  explained.  Hummel’s 
data  pertain  to  a  sharp,  flat  wing  at  a  =  20.5  deg,  and  Rec  =  0.9  x  106  in  low-speed  flow. 
The  “salient  edge”  concept  clearly  should  be  applicable  in  that  case,  and  vortex  breakdown 
had  not  occurred.  The  data  in  Fig.  13  pertain  to  an  aircraft  at  =  0.6,  with  a  finite 
leading-edge  radius,  a  =  6  deg,  and  Rec  =  20  -  40  x  IQ6.  The  concept  of  a  salient  edge  is 
at  least  questionable  in  this  case,  and  one  would  have  to  suspect  that  the  higher  Reynolds 
number  brought  about  vortex  breakdown  and  the  loss  in  suction  shown  in  Fig.  13. 

A  finite  leading-edge  radius  also  characterizes  the  model  in  Fig.  15.  Even  though  the 
angle  of  attack  is  10  deg,  it  appears  that  the  flow  is  approaching  an  attached  condition  as  Re 
increases.  The  low  Mach  number  assures  that  shock-induced  separation  is  not  a  factor. 

In  the  context  of  this  report,  it  is  inferred  that  a  scale  effect  should  be  expected  if  the 
higher  Reynolds  number  that  one  wishes  to  simulate  is  sufficient  to  change  the  basic 
character  of  the  flow.  That  can  occur  if  shock  interaction  is  significantly  changed  or  if 
vortex  breakdown  is  precipitated  as  Reynolds  number  increases.  It  will  sometimes  be  feasible 
to  approximate  the  subsonic,  high-Re  limit  by  computing  the  inviscid,  attached-flow  values 
of  the  aerodynamic  coefficients. 


It  is  well  known  that  fuselages  or  missiles  at  high  angles  of  attack  may  experience  yawing 
moments  and  side  forces  created  by  asymmetrical  vortices  shed  from  the  inclined  body. 
Evidence  of  Reynolds  number  effects  has  been  presented,  cf.  Keener  et  al.  (Ref.  25)  or 
Jorgensen  (Ref.  26).  The  data  show  that  quite  unexpected  effects  are  likely  to  arise  from 
attempts  to  simulate  higher  Reynolds  numbers  by  using  boundary  layer  trips.  Figure  16, 
taken  from  Ref.  26,  shows  that  increasing  Reynolds  number  affected  side  forces  in  a 
dramatic  fashion.  However,  either  longitudinal  strips  of  roughness  or  grit  applied  to  the 
nose  of  a  body  of  revolution  have  been  found  to  diminish  side  forces  or  yawing  moments 
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when  they  are  located  so  as  to  cause  symmetrical  separation  and  vortex  production.  It  is 
obvious  that  simulation  of  higher  Reynolds  numbers  by  the  use  of  trips  almost  surely  would 
not  have  yielded  the  same  results  as  are  shown  in  Fig.  16,  and  it  is  clearly  dangerous  to  use 
trips  in  any  case  where  vortex  strength  and  location  may  be  directly  influenced  by  proximity 
of  the  trips.  Forces  and  moments  arising  from  asymmetric  vortex  patterns  often  will  have 
unpredictable,  random  signs  depending  on  minute  model  construction  flaws,  support 
misalignment,  or  free-stream  flow  asymmetries.  Testing  methods  and  data  interpretation  are 
especially  critical  when  this  type  of  flow  exists. 

Leading-Edge  Separation  on  3-D  Wing 
M  =  0.3  a  =  10  deg 
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Figure  15.  Variation  of  lift  coefficient  with  Reynolds  number 
where  leading-edge  separation  and  vortical  flow 
exist.  From  Waggoner  (Ref.  24) 
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Figure  16.  Effect  of  Reynolds  number  on  side-force  coefficient 
for  an  ogive-cylinder  body  of  revolution.  From 
Jorgensen  (Ref.  26),  credited  to  Keener  (Ref.  25) 
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Figure  16.  Concluded. 
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3.0  SIMULATION  REQUIREMENTS 

Earlier  investigations  of  techniques  for  coping  with  transonic  SBL-simulation  deficien¬ 
cies  have  led  to  various  recommendations  and  warnings.  In  general,  there  is  agreement  that, 

An  essential  requirement  of  model-scale  testing  is  to  reproduce  the  full-scale 
location  of  the  shock  wave.  (Pankhurst,  Ref.  27) 

To  that  end,  it  has  been  recommended  that  first  priority  be  given  to  making  the  boundary 
layer  turbulent  at  the  foot  of  the  shock  (if  that  full-scale  condition  is  expected).  Recalling 
that  Pearcey  et  al.  (Ref.  2)  defined  Model  A  and  B  flows  which  have  differing  simulation 
requirements,  we  note  that  the  principal  requirement  for  Model  A  flows  is  that  the  turbulent 
boundary  layer  condition  be  provided.  Transition  must  be  complete  upstream  of  the  shock, 
but  it  need  not  be  as  far  forward  as  it  is  in  full-scale.  This  will  allow  closer  approximation  to 
(6/ c) fuii-scaie  by  making  (x,/c)modei  >  (xt/c)fuii-scaie  even  though  Rec  modei  <  Rec  mu-scaie- 

The  dismissal  of  Reynolds  number-related  effects  in  the  cases  of  unseparated  flows  and 
Model  A  flows  with  bubble  separations  may  be  permissible  in  regard  to  shock  position,  but 
there  are  Reynolds  number  effects  on  skin  friction  and  wake  which  should  not  be  ignored 
(cf.  Refs.  7, 10,  and  28).  This  is  summarized  as  follows: 

As  long  as  there  is  no  separation,  the  behaviour  of  the  boundary  layer  in  the 
region  of  interaction  does  not  appear  to  have  an  appreciable,  direct  effect  on  the 
position  of  the  shock  wave.  Also,  the  difference  between  the  actual  pressure 
distribution  and  one  with  a  discontinuous  pressure  rise  at  the  shock  position  has 
little  effect  on  the  lift,  drag,  and  pitching  moment  of  the  aerofoil.  However, 
boundary-layer  behaviour  through  the  interaction  does  have  an  effect  on  the 
subsequent  development  of  the  boundary  layer  and  wake  and  hence,  via  the 
weak,  overall  interaction,  it  does  have  some  influence  on  the  flow  field  at  large. 
(Green,  Ref.  7) 

If  separation  is  not  present,  the  shapes  of  the  velocity  profiles  in  the  boundary  layers 
approaching  the  trailing  edge  on  the  upper  and  lower  surfaces  are  the  key  features  to 
simulate,  according  to  Green.  In  an  investigation  reported  by  Blackwell  (Ref.  29)  where  a 
“small  amount”  of  shock-induced  separation  was  said  to  exist,  substantial  Reynolds 
number  effects  were  found,  and  the  recommendation  was  to  simulate  the  boundary  layer 
characteristics  at  the  airfoil  trailing  edge.  Matching  the  shape  factor,  H,  was  specifically 
recommended. 
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If  Model  B  flows  are  involved,  a  good  summary  of  the  situation  is  the  following. 

The  prime  requirement  for  correct  simulation  thus  remains  a  boundary  layer 
that  is  turbulent  at  the  point  of  interaction  with  the  shock.  It  is  now  also  clear, 
however,  that  there  is  a  further  requirement  of  almost  equal  importance  and 
generality:  this  turbulent  layer  must  be  of  a  thickness  that  is  not  so  magnified  in 
relation  to  full  scale  (and  must  have  a  profile  that  is  not  so  distorted)  that  after 
interaction  with  the  shock  it  will  provoke  a  rear  separation  that  would  not  be 
present  at  full  scale,  or  will  increase  significantly  the  severity  of  such  rear 
separation  as  might  be  present  at  full  scale.  (Pearcey  et  al.,  Ref.  2) 

The  latter  quotation  is  only  a  portion  of  the  “Concluding  Remarks”  in  Ref.  2,  which  should 
be  read  in  its  entirety  for  a  valuable  lesson  on  simulation  problems  of  Model  A  and  B  flows. 

Leaving  aside  the  question  of  how  one  is  to  predict  where  the  critical  airfoil  station  may 
be  in  a  given  case,  it  still  may  not  be  possible  to  bring  about  the  desired  turbulent  boundary 
layer  thickness  by  using  trips.  Tripping  can  be  supplemented  by  suction  or  cooling  to  thin 
the  tripped  layer,  but  the  expense  is  obvious  and  there  would  remain  a  problem  of 
accommodating  varying  test  conditions  without  an  elaborate  boundary  layer  control  system 
such  as  Green  and  Bore  have  discussed  (Refs.  7  and  30,  respectively). 

The  situation  faced  by  experimenters  who  cannot  meet  all  of  the  applicable  simulation 
requirements  necessitates  that  extra  steps  in  testing  and  computation  be  taken  if  full-scale 
results  are  to  be  predicted.  The  procedure  obviously  is  fraught  with  uncertainty,  and  there 
may  be  angles  of  attack  or  other  conditions  where  no  approach  is  going  to  be  successful. 
Figure  17  shows  the  Reynolds  number  gap  that  may  have  to  be  bridged  by  extrapolation. 

Basically,  if  subscale  measurements  can  be  made  at  conditions  where  reasonably  well- 
behaved  variations  with  increasing  Reynolds  number  can  be  relied  on,  then  theoretical  and 
computational  means  can  be  used  to  perform  the  required  extrapolation  to  higher  Reynolds 
numbers.  To  carry  this  out,  fuil-scale  position  of  the  shock  wave(s)  will  have  to  be  estimated, 
and  boundary  layer  tripping  usually  will  have  to  be  used  to  approximate  the  desired 
boundary  layer  conditions  at  the  shock.  Transport  aircraft  wings  will  be  affected  mainly  by 
Model  A  and  B  SBL-interaction  effects,  while  the  lower  aspect  ratios,  higher  angle-of-attack 
operations,  thinner  airfoils,  and  greater  sweep  of  fighters  and  missiles  often  will  make 
vortex-shock  interactions  important  in  simulation  of  those  vehicles.  Missiles  and  airplane 
forebodies  at  angles  of  attack  also  experience  flow  separation  and  vortex  formation,  with 
possible  scaling  problems  in  their  testing  (cf.  Fig.  16). 


29 


AEDC-TR-84-23 


In  summary,  the  currently  recognized  requirements  for  simulating  transonic,  viscous 
flows  are  as  follows  (assuming  equality  of  Mach  numbers  and  accurately  scaled  models): 

1.  Create  full-scale  shock  position,  i.e.,  (x/c)moije|  sh  =  (x/c)fui|.scaie  sh 

2.  Have  turbulent  boundary  layer  at  foot  of  shock  if  that  is  the  full-scale  condition 

3.  Have  (5  / c)model  ~  full-scale 

4.  Have  Hm0|jc|  *  H?uh  -scale 

5.  Create  full-scale  vortex  relative  position  and  strength  if  that  is  a  feature  of  the 
flow 

where  S*  =  boundary  layer  displacement  thickness;  c  =  wing  chord  or  other  appropriate 
length;  and  H  =  boundary  layer  shape  factor. 


at  shock 
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In  this  case,  the  “approximately  equal”  signs  must  be  regarded  as  meaning  “sufficiently 
close  to.”  In  typical  aerodynamic  testing,  these  requirements  are  seldom  all  met,  and  the 
foregoing  list  does  not  recognize  the  evidence  that  the  unit  Reynolds  number  affects  SBL 
interaction.  Of  course,  equality  of  model  and  full-scale  unit  Reynolds  numbers  in  scale- 
model  testing  is  incompatible  with  (1)  and  (3).  Testing  over  a  range  of  U/»  and  varying 
model  sizes  would  seem  the  only  recourse.  Number  (3)  becomes  virtually  impossible  when 
Rec  mwjel  <  <  Rec  fuLi-scale  and  boundary  layer  tripping  is  done,  unless  some  auxiliary  means 
of  thinning  the  boundary  layer  is  implemented.  Number  (4)  seems  within  reach  in  many 
cases.  Number  (5)  refers  to  a  complex  flow  interaction  which  may  be  influenced  by  very 
small  model  or  free-stream  asymmetries  such  that  the  flow  sometimes  seems  to  have  a 
random  character.  Requirements  numbered  (2),  (3),  and  (4)  have  generally  been  considered 
the  keys  for  satisfying  requirement  (1)  and  the  type  and  extent  of  separation  related  to  the 
SBL  interaction  on  wings  at  low  angles  of  attack.  It  is  desirable  to  simulate  full-scale 
boundary  layer  characteristics  at  the  trailing  edge,  as  well  as  at  the  shock,  but  one  may  not 
be  able  to  do  both.  Meeting  the  listed  conditions  at  the  highest  available  (U/iO*,  may  result 
in  the  closest  approach. 


It  is  important  to  be  aware  that  significantly  different  flow  conditions  can  exist  within 
the  general  class  of  transonic  flows.  Small  changes  in  Mach  number  can  strongly  influence 
viscous  interaction;  boundary  layer  transition  is  differently  affected  by  sweep  angle  at 
different  Mach  numbers;  roughness  rapidly  declines  in  effectiveness  as  Mach  number 
increases  in  the  region  of  Mk  «  1;  etc.  Therefore,  the  experimenter  seeking  to  establish 
dynamic  similarity  has  to  be  careful  of  generalizations  and  be  mindful  of  the  particular  case 
at  hand. 

Nothing  has  been  said  about  the  situation  that  will  arise  when  vortices  and/ or  wakes 
from  canard  surfaces  interfere  with  wings.  Whether  unusual  Reynolds  number  effects  arise 
in  such  cases  is  not  clear;  no  data  have  been  noted.  However,  it  is  not  difficult  to  visualize 
qualitative  Reynolds  number  effects.  Considering  that  configurations  with  gaps  are 
sometimes  encountered,  it  should  be  mentioned  that  equal  (gap  width/5*)  seems  more 
important  than  geometrically  scaled  gap  width  in  subscale  testing. 

4.0  INFLUENCING  BOUNDARY  LAYERS  TO  SIMULATE 
HIGHER  REYNOLDS  NUMBERS 

4.1  BOUNDARY  LAYER  TRANSITION  AND  THICKNESS 

There  have  been  recent  and  comprehensive  reviews  of  boundary  layer  transition  and  the 
nature  of  disturbances  created  by  single  roughness  elements.  For  that  information,  readers 
may  turn  to  Refs.  31,  32,  and  33.  The  specific  concern  with  transition  that  is  examined  here 
is  its  manipulation  so  as  to  compensate  for  the  deficiency  in  Reynolds  number.  The  major 
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object  is  to  bring  about  the  shock  location  and  boundary  layer  separation  that  duplicate 
what  would  occur  in  full-scale  free  flight.  Often  that  requires  moving  transition  forward  to  a 
station  where  it  would  occur  naturally  if  the  Reynolds  number  were  greater.  That  relocation 
may  be  done  by  using  boundary  layer  trips,  but  trips  can  cause  excessive  thickening  of  the 
boundary  layer  and  prevent  reproducing  full-scale  separation  patterns.  Trips  may  also  cause 
extraneous  drag  when  they  are  too  large.  Thus,  tripping  the  boundary  layer  at  the  desired 
location  is  not  the  only  requirement;  it  must  be  done  with  minimal  loss  of  energy  and 
extraneous  drag. 

Peterson  (Ref.  34)  measured  boundary  layer  profiles  relatively  far  downstream  of  trips 
(x/k  =  430  to  4560),  with  results  such  as  are  shown  in  Table  1.  Narrow-band  grit  roughness 
is  represented  in  these  data,  which  suggest  that  thickening  of  the  boundary  layer  is 
insignificant  when  k/S*  <!  1/2  and  x  >  >  k.  It  is  also  encouraging  to  note  that  the  shape 
factor  is  unaffected,  even  when  k  >  5k-  Subscript  1  means  no  trip;  2  means  with  trip. 


Table  1.  Data  on  Supersonic  Boundary  Layer  Downstream 
of  3-D  Trips.  From  Peterson  (Ref.  34) 


Me  *  3.1,  =  0.015  cm,  xk  =  0.64  cm,  x  = 

21.6  cm 

k/5k 

62/61 

55/5  T 

02/0| 

H2/Hi 

0 

1.000 

1.000 

1.000 

1.000 

0.307 

1.013 

0.980 

1.000 

0.983 

0.527 

1.026 

0.970 

0.980 

0.983 

0.607 

1.126 

1.157 

1.155 

1.000 

1.033 

1.168 

1.147 

1.150 

1.000 

1.293 

1.201 

1.265 

1.263 

1.000 

1.907 

1.226 

1.265 

1.269 

1.000 

1.987 

1.278 

1.304 

1.314 

0.983 

3.273 

1.291 

1.147 

1.184 

0.966 

P”  Xk  ^ 


However,  it  must  be  emphasized  that  these  data  represent  3-D  trips,  Mach  3  flow,  x  >  > 
xk  and  x  >  >  k.  Viewing  the  low-speed  data  of  Klebanoff  and  Tidstrom  (Ref.  35)  reveals 
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that  the  boundary  layer  at  stations  only  a  short  distance  downstream  of  the  roughness  is 
affected  by  even  very  small  trips.  Figure  18  from  Ref.  35  makes  this  point.  The  boundary 
layer  thicknesses  of  the  undisturbed  laminar  layer  at  the  roughness  station,  xk,  are  not  given 
in  Ref.  35,  but  calculated  values  for  a  Blasius  layer  are  5k  =  0.843  cm,  6k  =  0.279  cm  (0.11 
in.),  and  0k  =  0.108  cm  (0.04  in.)  for  the  condition  corresponding  to  the  square  symbols. 
The  shape  factor,  H,  would  be  2.59.  Thus,  Fig.  18  shows  that  S’  immediately  downstream  of 
xk  is  first  greater  than  the  untripped  value;  then  it  drops  below  that  level,  and  it  appears  that 
it  begins  to  increase  when  x  >  0  (100  k).  Momentum  thickness,  0,  just  downstream  of  xk  is 
lower  than  its  undisturbed  level,  and  it  slowly  increases  toward  that  level  with  increasing  x. 
The  shape  factor,  H,  initially  higher  than  the  undisturbed  laminar  value,  has  attained  the 
undisturbed  value  when  (x  -  xk)  >  0  (100  k).  Most  transonic  testing  requires  that  xt  -  xk. 

The  more  general  studies  of  tripping  have  concerned  2-D  flows  with  no  pressure 
gradients;  thus,  the  results  may  not  be  directly  applicable  to  cases  encountered  in  testing 
many  aircraft  models.  Reexamining  Fig.  7  arouses  suspicion  that  tripping  criteria  for  swept 
3-D  wings  differ  from  2-D  cases.  A  few  of  the  possible  factors  in  tripping  3-D  boundary 
layers  are  briefly  examined. 

Regarding  swept  wings,  there  is  some  reason  to  hope  that  2-D  methods  for  sizing  trips 
can  be  adopted  successfully  because  of  the  independence  principle.  Based  on  the  analysis  of 
laminar,  incompressible  flows  over  yawed  2-D  bodies,  Prandtl  (Ref.  36),  Struminsky  (Ref, 
37),  Jones  (Ref.  38),  and  Sears  (Ref.  39)  have  concluded  that  the  streamwise  flow  is 
independent  of  the  spanwise  flow  in  such  cases.  There  is  experimental  evidence  that 
independence  also  prevails  for  the  turbulent  boundary  layer  (Ref.  40),  but  others  have 
reached  the  opposite  conclusion  (Ref.  41).  However,  it  is  laminar  flow  that  has  to  be 
tripped,  and  the  independence  principle  offers  justification  for  using  the  same  method  of 
trip  sizing  on  both  swept  and  unswept  wings.  It  also  offers  encouragement  that  methods 
such  as  Refs.  42,  43,  44,  45,  and  46  can  be  used  if  confined  to  the  case  where  transition  is 
wanted  "near”  the  trip,  i.e.,  xt  -  xk.  The  latter  statements  are  made  even  though  the 
independence  principle  is  not  valid  in  boundary  layer  stability  theory  (Personal 
communication  from  Prof.  Eli  Reshotko  of  Case-Western  Reserve  University,  November 
1983).  It  is  possible  that  stability  considerations  do  not  apply  when  a  strong  overriding 
disturbance  such  as  trips  brings  about  transition. 

Evans  (Ref.  47)  examined  the  effectiveness  of  trips  on  several  swept  wings  and  found  his 
data  did  not  fit  the  predictions  for  2-D,  zero-pressure-gradient  flows  of  Refs.  42,  43, 44,  and 
45.  Specifically,  the  roughness  was  less  effective  than  predicted  by  the  methods  of  Refs.  43, 
44,  and  45,  but  more  effective  than  predicted  following  Ref.  42.  There  are  not  enough 
details  given  in  Ref.  47  to  permit  a  more  critical  examination  of  the  results.  For  example,  the 
pressure  distributions  on  the  wings  are  not  given.  Nash  and  Bradshaw  (Ref.  48)  have  shown 
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that  drag  of  excrescences  on  bodies  such  as  airfoils  can  be  much  greater  than  the  same 
roughness  on  a  flat  plate.  This  suggests  that  the  tripping  power  of  roughness  may  also  be 
greater  on  airfoils  where  local  velocities  exceed  free  stream. 


•5246 

x  -  x^.  In. 


Figure  18.  Distribution  of  boundary  layer  parameters  downstream 
of  0.06<Mn.-diam  roughness  elements  at  xk  =  2  ft. 
From  Klebanoff  and  Tidstrom  (Ref.  35) 
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The  Nash-Bradshaw  analysis  presents  an  insight  or  framework  for  understanding  the 
manner  in  which  a  perturbation  in  boundary  layer  momentum  thickness  adds  to  profile 
drag.  They  considered  only  boundary  layers  that  are  already  turbulent  at  the  roughness 
element.  However,  if  it  is  assumed  that  transition  occurs  in  a  short  distance  downstream  of 
the  roughness,  the  excrescence-drag  magnification  factor  due  to  a  non-flat-plate  flow  is 
applicable.  To  estimate  excrescence  drag,  it  is  necessary  to  have  a  value  for  the  Ad  added  by 
the  roughness  (cf.  Fig.  18). 

To  gain  some  appreciation  of  the  relative  sizes  of  boundary  layer  trips  that  may  be 
necessary  in  a  transonic  flow,  results  of  a  simple  example  are  tabulated  below.  Because  of  its 
convenience,  plus  the  present  concern  with  transition  near  the  trips,  the  method  of  Van 
Driest  and  Blumer  (Ref.  45)  is  followed.  The  calculated  trip  height,  k,  is  said  to  correspond 
to  the  minimum  transition  Reynolds  number.  However,  the  k  is  not  usually  large  enough  to 
bring  about  transition  at  the  trip  (xt  «  xk).  This  distinction  must  be  made  because  the 
decrease  of  xt  as  k  or  U/ v  increases  usually  slows  as  Xt  — ■  xk  and  a  minimum  Ret  for  a  given  k 
results.  Therefore,  it  is  better  to  view  the  k  calculated  by  the  method  of  Van  Driest  and 
Blumer  or  Braslow  et  al.  as  a  trip  height  that  gains  the  most  forward  movement  of  xt  that  can 
be  accomplished  “easily,”  i.e.,  without  encountering  the  stiffening  resistance  that  arises 
when  there  is  an  effort  to  move  xt  too  near  to  xk.  The  procedure  given  by  Potter  and 
Whitfield  (Ref.  42)  has  the  advantages  of  predicting  the  movement  of  xt  for  any  xk  <  xt  < 
xto,  up  to  the  case  xt  »  xk,  and  the  use  of  the  untripped  transition  station,  xt0  in  the 
correlation  removes  the  influence  of  “wind-tunnel-peculiar”  conditions  which  affect 
transition.  However,  a  more  tedious  computation  is  necessary  when  k  <  5k,  and  the  value  of 
Reynolds  number  of  transition  without  roughness  is  required. 

Because  the  method  of  Ref.  42  attempts  to  encompass  the  case  xt  *=  xk,  it  often  has  been 
used  to  determine  a  value  of  k  for  that  case,  when  the  user  of  the  method  would  have  had 
better  results  by  only  taking  k  large  enough  to  bring  about  the  decrease  of  xt  to  the  ‘  ‘near 
xk”  station  such  as  the  “effective  x:”  of  Ref.  45.  This  would  follow  from  limiting  the 
tripping  parameter  to  Re£/e.  =  0.8  -  0.85  rather  than  1.0  in  the  method  of  Ref.  42. 

To  furnish  an  example  of  the  sizes  of  trips  that  may  be  required  in  a  transonic  wind 
tunnel  case,  the  simpler,  though  less  general,  method  of  Ref.  45  is  applied  to  representative 
(hypothetical),  2-D,  flat-plate  flow  conditions  specified  in  Table  2.  The  results  given  there 
are  presented  as  a  ready  reference  to  the  range  of  k/6k  encountered  under  typical  transonic 
wind  tunnel  conditions. 

There  are  circumstances  which  will  alter  the  required  sizes  of  trips  or  their  effectiveness. 
Among  those  is  the  local  flow  field  such  as  already  mentioned  in  regard  to  the  effect  of 
pressure  gradients.  In  addition,  it  should  be  noted  that  Ref.  42  shows  that  when  Mk  <  1,  a 
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2-D  roughness  element  is  more  effective  than  the  3-D  grit  roughness  for  which  Ref.  45  and 
Table  2.  apply.  This  enables  a  much  lower  k/5k  in  subsonic  flow  if  a  2-D  element  meets  other 
requirements.  Secondly,  sweep  apparently  decreases  transition  Reynolds  numbers  as 
demonstrated,  e.g.,  by  Jillie  and  Hopkins  (Ref.  49).  They  found  that  free-stream  Reynolds 
number  of  transition  was  decreased  more  than  could  be  accounted  for  by  the  shock  loss  at  a 
blunt  leading  edge.  Their  data  for  =  0.27,  where  no  leading-edge  shock  could  exist,  are 
shown  in  Fig.  19.  Figure  20  shows  the  effect  of  sweep  when  combined  with  the  effect  of 
leading-edge  bluntness  and  normal  shock  formation  at  Mm  =  2.  Greater  reductions  in 
Reynolds  number  were  found  in  the  case  of  supersonic  flows,  where  leading-edge  shock  loss 
is  decreased  as  sweep  increases.  (The  dashed  curve  in  Fig.  20  represents  the  computed  effect 
of  a  normal  shock  wave  on  local  Reynolds  numbers.)  This  situation  will  make  for  easier 
boundary  layer  tripping;  i.e.,  less  trip  height  will  suffice  on  swept  wings.  The  marked  drop  in 
Reynolds  number  ratio  and  occurrence  of  separation  when  the  leading-edge  flow  component 
became  subsonic  in  Fig.  20  should  be  noted. 


Table  2.  Example  Trip-Size  Estimates  (Method  of  Ref.  45) 


Assumed  Conditions,  Me  =  1.00,  (\J/v)c  —  100,000/cm, 

Tw  -  Taw.  3-D  Roughness  Elements 

xk  cm 

Reck 

k,  cm 

k/5k 

xt,  cm 
(Predicted) 

1 

920 

0.009 

0.44 

1.9 

2 

1100 

0.011 

0.37 

2.9 

4 

1300 

0.013 

0.31 

4.9 

8 

1550 

0.016 

0.27 

8.9 

16 

1840 

0.018 

0.22 

16.9 

The  possibility  of  tripping  in  many  cases  without  exceeding,  say,  k/Sk  *  1/2  suggests 
that  excessive  boundary  layer  thickening  attributable  to  trips  should  not  be  as  severe  a 
problem  as  widely  assumed,  provided  one  is  concerned  with  stations  where  x  >  >  k.  Of 
course,  without  more  data  such  as  those  in  Table  1 ,  it  may  be  wrong  to  conclude  that  k/5k  * 
1/2  is  a  safe  limit  even  when  x  >  >  k.  Mach  number,  pressure  gradient,  type  of  trip,  lateral 
spacing  of  trip  elements,  and  existence  or  nonexistence  of  spanwise  flow  may  affect  the 
distortion  of  the  boundary  layer  downstream  of  the  trip  and  the  incremental  drag  due  to  the 
trip.  That  part  of  the  added  drag  caused  by  the  trip  rather  than  the  tripping,  i.e. ,  excrescence 
drag,  has  been  measured  by  various  investigators.  Braslow  (Ref.  43)  and  Winter  and  Gaudet 
(Ref.  50)  have  reported  such  data.  The  latter  also  gave  an  analysis  and  method  of  predicting 
excrescence  drag.  Braslow  reported  cases  where  k  =  did  not  cause  added  drag,  6k  being 
laminar  boundary  layer  thickness.  However,  that  is  not  a  general  criterion.  Trip-added  drag 
apparently  is  magnified  when  local  velocities  exceed  free-stream  values,  cf.  Nash  and 
Bradshaw  (Ref.  48). 
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Figure  19.  Variation  of  the  normalized  frcc-stream  transition  Reynolds 
number  with  sweep  at  subsonic  Mach  number.  From  Jillie 
and  Hopkins  (Ref.  49) 


There  is  another  aspect  of  the  boundary  layer  problem  that  arises  even  when  the  trip  adds 
little  to  the  thickness.  A  low  Reynolds  number  boundary  layer  that  is  laminar  over  a  distance 
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of  approximately  xk  and  turbulent  thereafter  can  equal  the  relative  thickness  (£/c)  of  a  high 
Reynolds  number  layer  that  is  naturally  turbulent  from  xt  *  0  only  under  certain  conditions. 
Matching  of  relative  thicknesses,  5*/c,  0/c,  at  equal  values  of  (x/c)s,  shock  location,  between 
a  wind  tunnel  case  and  a  higher-Reynolds-number  flight  case  depends  on  the  satisfaction  of 
several  flow  and  geometric  relations  which  is  not  always  feasible. 


Sweep,  A  ,  deq 
*>  2. 00 


Figure  20.  Variation  of  the  normalized  free-stream  transition  Reynolds 
number  with  sweep  at  supersonic  Mach  number.  From 
Jillie  and  Hopkins  (Ref.  49) 
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The  trade-off  of  the  effects  of  local  Mach  number  and  local  Reynolds  number  at  the 
boundary  layer  edge  should  be  kept  in  mind  when  choosing  trip  location.  The  rapid  increase 
in  boundary  layer  stability  through  the  transonic  Mach  number  range  causes  a 
corresponding  rate  of  increase  in  trip  sizes  if  unit  Reynolds  number  is  not  also  increasing 
enough  to  compensate  for  the  Mach  number  effect,  cf.  Ref.  42. 

4.2  VORTEX  GENERATORS 

The  use  of  vortex  generators  to  inhibit  boundary  layer  separation  is  an  established 
practice.  Figure  21 ,  from  Cahill  (Ref.  51),  is  an  interesting  example.  With  transition  fixed  at 
x,/c  =  0.10  and  a  «  0  deg,  shock  location  and  trailing-edge  pressure  at  Rec  =  2.8  x  106 
were  nearly  the  same  as  found  for  Rec  =  8.5  x  106  if  vortex  generators  were  located 
upstream  of  the  shock  (diamond  symbols  in  Fig.  21).  When  the  shock  moved  ahead  of  the 
vortex  generators,  at  a  =  1  to  2  deg,  the  trailing-edge  pressure  and  shock  location 
dramatically  shifted  to  the  levels  corresponding  to  lower  Reynolds  numbers.  That 
underscores  the  problem  with  either  trips  or  vortex  generators  at  fixed  locations.  Transition, 
shock  waves,  and  separation  move  about  as  Mach  number  and/ or  angle-of-attack  change. 
Thus,  no  single  trip  or  vortex  generator  position  is  likely  to  satisfy  all  test  requirements.  This 
problem  was  encountered  long  ago,  and  for  the  trips,  some  degree  of  flexibility  sometimes 
has  been  gained  by  doing  the  tripping  with  gas  jets  that  could  be  varied  in  stength. 
Additional  flexibility  also  may  be  achieved  by  arranging  to  exhaust  the  jets  at  varying 
stream-wise  locations,  but  that  obviously  requires  a  more  complex  model. 

Under  some  circumstances,  it  may  be  feasible  to  combine  the  vortex  generation  and  trip 
function,  i.e.,  use  the  generators  as  trips.  Some  trade-off  in  effectiveness  in  fulfilling  each 
function  can  be  expected.  For  example,  the  higher  velocities  near  the  surface  in  a  turbulent 
boundary  layer  probably  enhance  vortex  generator  performance,  but  the  thicker  turbulent 
boundary  layer  may  offset  that  advantage. 

The  design  of  vortex  generators  has  been  discussed  by  Powers  (Ref.  52)  who  has  given 
criteria  for  generator  angle  of  attack,  spacing,  aspect  ratio,  and  height.  The  function  of 
vortex  generators  is  to  cause  the  transport  of  higher  speed  air  from  the  outer  part  of  a 
boundary  layer  into  the  lower  speed  region  nearer  to  the  solid  surface.  When  adjacent 
generators  are  canted  so  that  counter-rotating  vortices  are  formed,  the  boundary  layer 
flowing  between  the  pair  can  be  made  much  thinner.  This  suggests  that  vortex  generators 
may  be  used  to  simulate  a  boundary  layer  at  a  higher  Reynolds  number  if  the  added  drag  of 
the  generators  is  acceptable.  Considering  that  a  generator  need  be  no  higher  than  6,  will  be  at 
an  incidence  of  approximately  16  deg  with  the  local  flow,  will  have  an  aspect  ratio  near  1, 
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Figure  21.  Effect  of  vortex  generators  on  shock  location  and 
pressure  recovery.  From  Cahill  (Ref.  51) 
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and  be  spaced  about  3.5  generator  heights  from  its  neighbors,  it  would  seem  that  the  extra 
“trip”  drag  would  be  no  worse  than  that  caused  by  rows  of  grit  or  other  conventional  trips. 
Unfortunately,  there  seems  to  be  no  data  on  the  development  of  transition  downstream  of 
vortex  generators,  so  their  performance  as  trips  cannot  be  compared  to  the  usual  2-  or  3-D 
trips.  Aside  from  the  added  work  in  making  such  devices,  another  disadvantage  is  the 
spanwise  nonuniformity  because  of  the  wake-like  flow  that  will  exist  upstream  of  the  station 
where  merging  occurs. 

4.3  TREE-STREAM  TURBULENCE 

Not  only  is  free-stream  turbulence  a  factor  in  determining  boundary  layer-transition 
Reynolds  number,  it  is  also  known  to  affect  boundary  layer  thickness  and  detailed 
characteristics  within  the  layer.  Raghunathan  and  McAdam  (Ref.  53)  have  presented  data 
for  turbulent  boundary  layer  flow  over  a  half-biconvex  airfoil  mounted  on  the  wall  of  a 
transonic  wind  tunnel.  Ratios  of  total  and  momentum  thicknesses  and  shape  factors  across 
the  shock  are  given.  Increasing  free-stream  turbulence  causes  62/^1  and  H2/Hi  to  decrease, 
while  62/6i  increases  rather  slowly.  Separated  flow  region  at  the  trailing  edge  was  decreased 
as  free-stream  turbulence  increased.  In  Ref.  54  the  same  authors  presented  data  for  the  same 
experimental  arrangement  which  show  shock  position  shifted  downstream  as  free-stream 
turbulence  was  increased. 

Green  (Ref.  55)  reviewed  earlier  data  on  flat-plate  flows  that  predicted  that  increased 
stream  turbulence  (at  constant  Re*)  will  alter  skin-friction  coefficients  and  shape  factors 
according  to  the  relations,  Acf/cf  =  4.8  u'/uE,  and  AH/H  =  -(2.4  -  0.25H)u'/u*.  Green 
concluded  that  the  “effective”  Reynolds  number  increases  rapidly  with  u'/ue  and  that  the 
sensitivity  is  greater  at  higher  Reynolds  number.  In  a  typical  example,  for  an  airfoil,  he 
found  that  a  1-percent  increase  of  turbulence  could  have  the  same  effect  as  increasing  the 
Reynolds  number  by  60  to  70  percent. 

These  findings,  not  surprisingly,  have  inspired  some  to  suggest  that  elevating  free-stream 
turbulence  is  a  way  to  gain  the  effects  of  higher  Reynolds  number.  However,  Green 
cautioned  that  there  were  unanswered  questions  about  the  validity  of  that  technique.  He 
pointed  out  that  the  role  of  turbulence  scale  as  well  as  intensity  would  have  to  be 
investigated;  the  response  of  the  boundary  layer  on  the  model  may  vary  depending  on  its 
scale;  and  acoustic  disturbances  may  not  affect  a  boundary  layer  in  the  same  way  as  vortical 
disturbances  coming  from  upstream.  Other  objectional  features  of  this  technique  that  can  be 
foreseen  are  the  need  to  calibrate  the  turbulence  and  its  equivalent  Reynolds  number  and  the 
detailed  relationship  between  turbulence  and  its  influence  on  model  measurements  even 
when  transition  is  not  a  factor  (cf.  Durbin  and  Hunt,  Ref.  56). 
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In  particular,  the  possibility  that  intensity  and  scale  of  turbulence  may  each  affect 
boundary  layer  flow  in  an  interrelated  fashion  cannot  be  ruled  out.  This  would  make  any 
effort  to  control  “effective”  Reynolds  number  by  manipulating  stream  turbulence  very 
risky.  It  opens  another  subject  of  research.  Of  course,  one  cannot  say  that  turbulence  scale 
effects  are  not  involved  in  conventional  boundary  layer  tripping  either.  Scale  of  trip- 
produced  turbulence  probably  is  related  to  trip  size,  local  Reynolds  number,  and  distance 
downstream  of  the  trip.  It  is  possible  that  scale  effect,  if  it  is  significant,  may  be  accounted 
for  by  the  usual  practice  of  evaluating  trip  effectiveness  in  terms  of  a  trip  height-Reynolds 
number  parameter. 


5.0  WIND  TUNNEL  TESTING  PRACTICES 

A  report  which  was  written  by  Haines  et  al.  (Ref.  57)  over  twenty-six  years  ago  presents  a 
rather  thorough  discussion  of  various  problems  encountered  in  fixing  transition  on  wings  in 
transonic  wind  tunnel  tests.  Most  of  the  noted  difficulties  remain  to  this  day.  A  few  of  the 
more  recent  reports  in  the  same  area  are  Refs.  11,  58,  and  59.  The  latter  also  covers  other 
transonic  wind  tunnel  testing  problems.  A  general  impression  left  by  all  such  discussions  is 
that  extra  effort  in  collecting  data  and  thoughtfulness  in  analysis  are  essential  if  the  best 
results  are  to  be  achieved  in  subscale  wind  tunnel  testing.  Routine  methods,  inflexible  test 
plans,  and  minimal  tunnel  occupancy  time  often  will  lead  to  later  surprises  when  full-scale 
data  become  available. 

Tripping  of  the  boundary  layer  is  not  necessarily  sufficient  to  assure  simulation  of  higher 
Reynolds  numbers,  the  problem  can  be  much  more  complex.  Different  unit  Reynolds 
numbers,  different  scales  (model  sizes),  different  trip  sizes  and  locations,  judicious  use  of 
vortex  generators,  guidance  by  computed  results,  measurement  of  boundary  layer 
characteristics  at  key  locations,  and  thoughtful  study  of  results  will,  in  general,  all  be 
required  in  order  to  minimize  the  error  when  the  test  results  are  extrapolated  to  full-scale 
conditions.  The  greater  uncertainity  generally  will  be  associated  with  3-D  flows  and  flows 
including  significant  shock  wave-vortex-transition  and  separation  interactions. 

If  wind  tunnel  testing  is  to  be  done  at  less  than  full-scale  Reynolds  number,  it  should  be 
guided  by  results  of  analytical  efforts  that  provide  estimates  of  full-scale  boundary  layer  and 
shock  wave  characteristics.  Otherwise,  the  testing  is  done  without  knowing  approximately 
the  conditions  that  will  indicate  how  well  full-scale  Reynolds  numbers  are  simulated.  The 
experimentation  will  include  determination  of  boundary  layer-transition  location  on  a 
smooth  model,  possibly  using  a  subliming  film  for  visual  location  of  transition  at  various 
free-stream  unit  Reynolds  numbers.  It  must  be  remembered  that  boundary  layer-transition 
Reynolds  number  usually  varies  with  unit  Reynolds  number,  and  other  properties  of  the  test- 
section  flow  of  transonic  wind  tunnels  also  often  vary  with  unit  Reynolds  number.  Thus, 
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Reynolds  number  effects  appearing  in  model  measurements  must  not  be  confused  with  those 
caused  by  free-stream  flow  changes  accompanying  changes  of  unit  Reynolds  number. 

Pressure  distribution,  or  at  least  trailing-edge  pressure,  and  boundary  layer  thicknesses 
at  the  trailing  edge  or  ahead  of  the  shock  may  also  be  measured  and  used  as  a  criterion  for 
scaling  viscous  effects.  Then,  it  is  usual  for  roughness  to  be  fixed  so  as  to  cause  transition  at 
location(s)  simulating  full-scale  boundary  layer  conditions  at  the  critical  location  on  the 
body.  Pressures  are  again  measured  to  determine  if  the  expected  full-scale  shock  position  has 
been  achieved. 

There  are  two  approaches  to  tripping,  which  may  be  classified  as  ‘’forward  tripping” 
and  ‘’aft  tripping.”  The  first  and  perhaps  more  traditional  forward  tripping  involves 
tripping  at  the  percent  of  chord  where  transition  is  predicted  under  full-scale  conditions; 
then  measurements  made  at  each  of  a  series  of  Reynolds  numbers  are  extrapolated  to  predict 
the  full-scale  results.  Aft  tripping  is  done  in  order  to  minimize  6Vc  at  the  shock  while  still 
assuring  a  turbulent  boundary  layer  there.  That  approach  seems  to  have  more  followers,  but 
it  cannot  always  be  used.  When  conditions  are  such  that  the  shock  or  separation  are  far 
forward,  tripping,  to  be  effective,  must  occur  still  farther  forward,  and  the  distinction 
between  forward  and  aft  tripping  may  vanish.  One  may  list  advantages  and  disadvantages 
for  both  approaches,  and  neither  will  necessarily  assure  that  trailing-edge  boundary  layer 
conditions,  which  affect  lift  and  drag,  are  closely  simulated.  Some  experimenters  will  choose 
to  use  each  method,  depending  on  the  particular  quantity  being  measured.  Any 
measurements  that  involve  significant  fore  and  aft  movement  of  separation  or  shock 
obviously  are  going  to  cause  problems  with  aft  tripping.  There  may  be  conditions  for  which 
"no  tripping”  is  the  preferred  approach.  Unfortunately,  it  may  be  impossible  to  forecast 
what  technique  is  best,  and  a  test  plan  should  allow  for  alternatives. 

If  there  is  concern  over  excrescence  drag,  the  drag  of  the  model  with  varying  trip  heights 
and  locations  may  be  measured.  If  possible,  schlieren  or  shadow  photography  may  be  used 
to  supplement  other  data.  Vortex  generators  also  may  be  used  to  delay  separation  and 
simulate  higher  Reynolds  number  conditions  on  the  model.  All  of  these  data  have  to  be 
collected  at  each  angle  of  attack.  The  pressure  distributions  may  be  used  in  the  method  of 
Ref.  15  to  estimate  shock  position  and  pressure  distribution  at  flight  Reynolds  numbers.  If 
full  pressure  distributions  are  not  taken,  the  trailing-edge  pressure  provides  the  indication  of 
which  trip  and/or  vortex  generator  configuration  has  achieved  the  highest  effective 
Reynolds  number. 

It  is  apparent  that  estimates  of  full-scale  boundary  layer  thickness  and  transition  location 
will  be  needed  if  the  effective  wind  tunnel  Reynolds  numbers  are  to  be  estimated.  That  is,  the 
tripped  boundary  layer  with  its  corresponding  characteristics  and  transition  location  has  to 
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be  related  to  an  untripped,  higher  Reynolds  number  condition.  Transition  location  alone 
could  be  used  as  the  scaling  criterion,  but  that  usually  is  inadequate.  More  detailed 
knowledge  of  the  tripped  layer  is  needed  because  thicknesses  and  shape  factors  offer  a  way 
to  judge  if  separation  is  also  properly  scaled. 

If  vortices  and  vortex-shock  interaction  are  features  of  the  flow,  it  will  be  necessary  to  try 
to  assure  that  the  line  of  separation  is  the  same  on  the  model  as  it  would  be  at  flight  Reynolds 
numbers.  This  may  require  that  leading-edge  roughness  be  used,  and  the  minimization  of 
boundary  layer  thickening  due  to  the  tripping  is  critical.  Various  computational  estimates  of 
separation  location  may  be  used  to  assess  Reynolds  number  effects.  The  interaction  of 
vortices  with  adjacent  surfaces  will  make  it  important  to  achieve  the  same  vortex  pattern  as 
will  exist  in  flight  (Ref.  59).  In  all  transonic  testing  the  criticality  of  Mach  number  as  a 
parameter  cannot  be  overemphasized.  That  is  also  a  troublesome  factor  when  attempting  to 
compare  flight  and  wind  tunnel  data  because  a  close  match  in  Mach  numbers  is  required 
when  near  Mach  one. 

It  is  important  that  the  type(s)  of  separation  be  identified  in  the  early  phase  of  testing  so 
that  the  essential  requirements  for  simulation  are  known.  This  information  is  also  needed  as 
a  check  on  the  applicability  of  Ref.  15  if  that  procedure  is  to  be  used. 

5.1  SCALING  FOR  AERODYNAMIC  COEFFICIENTS 

No  attempt  is  made  to  discuss  all  possible  scaling  procedures,  but  the  simpler  procedure 
appropriate  for  airfoils  or  wings  at  low  angles  of  attack  is  briefly  considered  as  an  example. 
We  concentrate  on  a  few  of  the  aerodynamic  qualities  that  have  to  be  scaled  when  transonic 
test  data  are  obtained  in  wind  tunnels  that  do  not  match  free-flight  Reynolds  numbers. 
Foremost  of  these  is  shock  location  on  wings  which  has  been  discussed  earlier  and  will  not  be 
reviewed  again.  However,  the  attainment  of  the  full-scale  shock  location  is  the  underlying 
basis  of  what  follows. 

Drag  coefficient  is  one  of  the  first  aerodynamic  coefficients  thought  of  in  the  context  of 
Reynolds  number  deficiency.  The  common  approach  is  to  acquire  data  over  the  available 
Reynolds  number  range  and  then  extrapolate  the  friction-drag  component  or  the  total  drag 
coefficient  to  the  desired  Reynolds  number  according  to  an  appropriate  friction-drag  law. 
“Form  factors”  are  used  as  corrections  to  flat-plate  friction  drag  to  account  for  the  shape  of 
the  airframe  components.  For  this  procedure  to  be  valid,  the  location  of  boundary  layer 
transition  must  either  be  fixed  or  its  variation  with  Reynolds  number  known,  and  no  further 
significant  change  in  shock  wave  or  separation  locations  can  occur.  Potential  errors  in  this 
procedure  are  easy  to  identify,  but  it  is  a  widely  used  technique,  and  the  only  avenue  toward 
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improvement  would  seem  to  be  work  on  details  of  the  basic  method.  For  example,  newer 
CFD  capabilities  will  improve  the  predicted  friction-drag  variation.  In  time,  CFD  will  enable 
the  entire  viscous  flow  field,  including  separation  to  be  scaled  more  precisely  on  a  routine 
basis. 

The  relatively  straightforward  approach  to  correcting  drag  coefficients  should  not  lead 
to  undue  optimism  about  extrapolation  procedures  in  general.  It  is  possible  for  movements 
of  transition  and/or  separation  to  cause  trouble.  A  sobering  example  is  presented  by 
Stanewsky  (Ref.  60).  His  data  show  a  reversal  of  trends  of  lift  coefficient  at  drag-rise  Mach 
number,  as  well  as  variation  in  that  Mach  number,  in  the  neighborhood  of  Rec  =  107  when 
free  transition  is  allowed.  With  transition  fixed  well  forward,  neither  Cl  nor  Ma  at  drag  rise 
agreed  with  those  values  or  even  their  trends  for  free  transition  at  the  maximum  Reynolds 
number  of  3  x  107. 

It  has  to  be  remembered  that  there  are  mutual  interactions  between  boundary  layer 
characteristics,  shock  locations,  and  separations,  and  these  are  varying  as  unit  Reynolds  and 
Mach  numbers  and  angle-of-attack  changes.  Thus,  there  can  be  unexpected  results  that  do 
not  follow  simple  extrapolations  of  lower  Reynolds  number  data.  One  must  keep  their 
attention  on  the  more  significant  parameters  that  influence  the  viscous-inviscid  interactions 
if  these  results  are  to  be  anticipated.  For  example,  if  one  wants  to  know  what  will  happen  if 
unit  Reynolds  number  is  increased,  they  have  to  be  mindful  of  the  effects  on  transition 
location  and  boundary  layer  characteristics,  the  influence  of  these  factors  upon  separation 
development  and  shock  location,  and  finally  the  net  effect  upon  pressure  distribution.  The 
interdependence  of  all  of  these  will  often  make  accurate  extrapolations  or  scaling 
impossible,  particularly  when  3-D  flow  is  involved. 

Stanewky’s  (Ref.  60)  research  has  shown  that  supercritical  airfoils  are  very  sensitive  to 
viscous  effects.  He  has  demonstrated  that  changes  in  unit  Reynolds  number  produced 
different  lift  coefficients  even  at  low  a,  and  that  reversals  in  movement  of  transition  location 
occurred  as  angle  of  attack  changed  because  of  accompanying  variations  of  boundary  layer 
thickness  and  effective  airfoil  shape.  Stanewsky’s  transonic  sensitivity  parameter  for  airfoils 
at  given  Ma  and  Rec  was  shown  in  Ref.  60  to  correlate  the  ratio  of  change  in  Cl  to  change  in 
transition  location  on  a  number  of  airfoils. 

One  approach  for  attempting  to  determine  aerodynamic  coefficients  for  an  airfoil  or 
wing  at  low  angle  of  attack  for  free-flight  Reynolds  number  may  be  based  on  the  attainment 
of  full-scale  shock  position  in  the  subscale,  low  Reynolds  number  tunnel  test.  See  earlier 
discussion  and  Fig.  22  for  the  justification  of  this  approach.  It  follows  from  conditions  (1) 
through  (4)  in  Section  3.0.  The  following  steps  would  be  involved,  for  given  Ma  and  a : 
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1 .  Estimate  transition  location,  (x/c)t,  for  no  trips  and  then  determine  trip  sizes  to 
bring  about  transition  at  a  series  of  lesser  (x/c)t.  The  intent  is  to  bring  about  a 
turbulent  boundary  layer  of  the  least  feasible  thickness  ahead  of  the  shock. 

2.  Measure  CL,  CD,  CM,  (x/c)sh  and  Cpte  at  each  (x/c),  at  the  maximum  chord 
Reynolds  number. 

3.  Plot  the  measured  data  as  in  Fig.  23.  If  the  rearmost  relative  shock  position, 
(x/c)sh  is  the  same  as  will  exist  at  full-scale,  the  coefficients  corresponding  to 
that  shock  position  should  be  approximately  the  ones  for  full-scale  Reynolds 
numbers. 

4.  Using  measured  Cpte,  apply  the  method  of  Ref.  15  to  further  check  on  shock 
position.  (The  limitation  on  that  method  relative  to  bubble  separations  has  to 
be  observed.)  Data  at  various  spanwise  stations  for  a  wing  model  and  for  other 
Reynolds  numbers  may  also  be  needed. 

5.  If  pressure  recovery  over  the  aft  part  of  the  chord  is  less  than  expected,  vortex 
generators  may  be  tested  at  several  chordwise  positions  upstream  of  the  shock 
to  learn  if  shock  position  is  moved  farther  aft  (see  Fig.  21  in  this  regard). 

6.  If  .vortices  are  a  feature  of  the  flow,  the  foregoing  procedure  should  be  carried 
out  both  for  vortices  burst  (broken  down)  and  unburst.  Tripping  may  not  be 
feasible  (see  earlier  discussion),  but  an  effort  has  to  be  made  to  at  least  bracket 
the  low  and  high  Reynolds  number  vortex  conditions  (see  Fig.  15). 

It  cannot  be  guaranteed  that  this  procedure  will  automatically  lead  to  the  correct  full- 
scale  aerodynamic  coefficients,  and  it  will  not  always  be  possible  to  approach  closely  enough 
to  conditions  (1)  through  (4)  of  Section  3.0.  Also  to  be  noted  is  the  need  for  more 
information  than  merely  the  low  Reynolds  number  wind  tunnel  data  to  carry  out  this 
procedure.  The  experimental  data  have  to  be  supplemented  by  CFD  and  correlations  such  as 
provided  in  Ref.  15. 

If  there  is  uncertainty  about  transition  location  and  no  preliminary  experiment  for  the 
purpose  of  locating  transition  is  conducted,  the  “en”  method  may  be  used.  A  recent  report 
by'  Wazzan  et  al.  (Ref,  62)  and  other  references  therein  may  be  referred  to  for  this 
computation.  However,  it  would  be  well  for  users  of  the  “en”  method  to  read  the  paper  by 
Mack  (Ref.  63).  Because  of  the  strong  destabilizing  influences  present  in  transonic  wind 
tunnels,  empirical  data  from  the  particular  tunnel  to  be  used  probably  will  be  relied  on  by 
many  experimenters.  Comparative  data  from  many  tunnels  have  been  reported  by 
Dougherty  and  Steinle  (Ref.  64)  and  Dougherty  and  Fisher  (Ref.  65).  A  recent  paper  by 
Harvey  (Ref.  66)  reviews  transonic  tunnel  and  flight  transition  results  in  relation  to  free- 
stream  disturbances. 
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Figure  23.  Idealized  graphical  trend  of  data  with  transition  location. 
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6.0  CONCLUDING  REMARKS 

When  this  subject  was  discussed  over  ten  years  ago,  all  of  the  parts  of  the  problem  were 
recognized.  However,  there  is  now  a  clearer  understanding  of  several  features  of  transonic, 
viscous-flow  simulation.  While  these  improvements  should  lead  to  better  degrees  of 
simulation,  they  do  not  make  the  experimenter’s  task  easier.  For  example,  it  has  been  shown 
that  when  unit  Reynolds  number  and  wetted  length  are  varied  separately,  different  effects 
on  separation  pressure  coefficient  are  produced  (Ref.  12).  This  type  of  effect  also  has  been 
demonstrated  in  the  upstream  influence  of  2-  and  3-D  compression  corners  (Ref.  13).  Thus, 
the  experimenter  cannot  be  satisfied  with  merely  matching  Reynolds  numbers  when 
investigating  shock  wave-boundary  layer  interactions. 

Melnik’s  viscous,  transonic  simulation  parameter,  Kt,  shows  that  Reynolds  number  (or 
cf)  affects  upstream  influence  of  shocks  even  though  the  boundary  layer  is  turbulent.  This 
imposes  an  extra  burden  on  scale-model  testing;  i.e.,  simply  achieving  a  turbulent  boundary 
layer  ahead  of  the  shock  is  not  sufficient  for  simulation  of  high  Reynolds  numbers. 

Another  aspect  of  the  simulation  problem  on  which  more  light  has  been  shed  is  the 
Reynolds  number  effect  on  separated  and  vortical  flows.  It  has  been  shown  that  changes  in 
aerodynamic  loads  because  of  vortices  may  persist  at  very  large  Reynolds  numbers  (Ref.  22), 
and  use  of  boundary  layer  trips  may  so  influence  vortex  formation  that  simulation  would  be 
entirely  lost  (see  Fig.  16).  This  serves  as  a  clear  warning  of  the  danger  of  assuming  that  all 
vortical  flows  are  immune  to  Reynolds  number  effects.  Shock  interaction  and  vortex 
bursting  may  be  influenced  by  large  changes  in  Reynolds  number,  and  the  possibility  of 
relatively  large  changes  in  aerodynamic  coefficients  should  not  be  ignored. 

Data  on  the  new  supercritical  airfoils  reveal  strong  sensitivity  to  viscous  effects. 
Stanewsky  (Ref.  60)  has  convincingly  documented  this. 

The  experimenter  is  compelled  to  recognize  several  general  features  of  viscous,  transonic- 
flow  simulation,  viz., 

•  The  technique  must  be  tailored  to  the  vehicle  configuration,  local  flow 
characteristics,  and  the  type  of  aerodynamic  data  sought. 

•  Tripping  the  boundary  layer  at  forward  stations  does  not,  in  general,  assure 
adequate  simulation  of  higher  Reynolds  numbers. 

•  Procedures  that  may  be  referred  to  as  “no  boundary  layer  tripping/’  “aft 
tripping,”  and  “forward  tripping”  may  each  have  its  place  in  a  test  program. 
Tripping  can  sometimes  conceal  Reynolds  number  effects. 
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•  Other  free-stream  flow  properties  of  transonic  wind  tunnels  often  vary  with 
changes  in  free-stream  unit  Reynolds  number,  thereby  creating  confusion 
between  “real”  and  “spurious”  viscous-flow  effects. 

•  Any  procedure  will  have  to  offer  guidance  and  alternatives;  it  cannot  be 
followed  blindly  without  incurring  risk. 

A  rational  subdivision  of  aerospace  vehicles  or  components  and  some  characterizations 
follow: 

1.  Transport  wings  —  low  angles  of  attack;  higher  Reynolds  numbers;  most 
extensively  studied;  great  accuracy  required  in  predicting  coefficients. 

2.  Delta  wings  —  vortical  flows  even  at  relatively  low  angles  of  attack;  may  be 
difficult  to  simulate  high  Reynolds  numbers  because  of  flow  developments  very 
near  leading  edges. 

3.  Missile  and  aircraft  forebodies  —  high  angles  of  attack;  vortical  flows;  difficult 
to  predict  results;  sometimes  large  effects  of  Reynolds  number. 

4.  “Complex”  configurations  —  a  category  of  complicated  flows  possibly 
involving  canards;  leading-edge  flaps;  relatively  large  inlet  and  nozzle-afterbody 
influences;  external  stores,  vectored  nozzles,  etc. 

Regarding  category  (1)  above,  the  current  questions  concern  whether  forward  tripping 
and  extension  of  wind  tunnel  data  to  predict  full-scale  aerodynamic  coefficients  by  simple 
trend  extrapolation  with  Reynolds  number  is  superior  to  aft  tripping  with  an  effort  to 
simulate  essentially  full-scale  flow  conditions  at  tunnel  Reynolds  number.  There  is  also 
concern  over  whether  matching  of  boundary  layer  properties  just  ahead  of  the  terminal 
shock  or  at  the  trailing  edge  is  best.  Actually,  all  of  these  have  their  places,  but  it  is  often 
unclear  which  plan  is  best  without  trying  all  of  them. 

For  category  (2),  it  will  sometimes  only  be  possible  to  bracket  the  low  and  high  Reynolds 
number  aerodynamic  coefficients.  A  particular  problem  arises  because  the  physical  space 
available  for  manipulating  the  boundary  layer  may  be  of  the  order  of  10  percent  of  the 
chord,  or  less,  at  the  leading  edge.  In  the  case  of  category  (3),  tripping  can  be  ineffective  or 
misleading.  The  many  different  and  complex  interactions  possible  in  category  (4)  make  a 
fixed  procedure  elusive.  None  is  suggested  herein. 

These  are  only  some  of  the  complexities  of  transonic,  viscous-flow  simulation  and 
scaling  that  have  been  brought  out  in  recent  years.  The  National  Transonic  Facility  (NTF)  at 
Langley  Research  Center  is  going  to  offer  much  greater  Reynolds  numbers  in  a  ground  test 
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facility,  so  that  an  increased  rate  of  learning  may  be  anticipated.  However,  much  testing  will 
continue  to  be  done  in  the  older,  conventional  tunnels  now  in  use.  This  review  stresses  the 
need  for  more  extensive,  systematic  experimental  procedures,  backed  up  by  CFD  and 
painstaking  analysis,  if  the  quality  of  information  derived  from  low  Reynolds  number  tests 
is  to  be  improved. 
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NOMENCLATURE 


Co  Drag  coefficient 

Cl  Lift  coefficient 

Cm  Pitching-moment  coefficient 

Cp  Pressure  coefficient 

Cps  Separation  pressure  coefficient 

Cp[e  Cp  at  trailing  edge 
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c  Chord  length  (c  =  average  chord) 

Cf  Local  skin-friction  coefficient 

H  Boundary  layer  shape  factor,  $70 

_  J 

H  (1/0  J  [Q/ee)(l  -  u/Ue)dy 

O 

K  Khan-Cahill  parameter,  Eq.  (6) 

Kt  Melnik’s  parameter  Eq.  (3) 

k  Height  of  roughness  element 

Lm  Upstream  influence  distance 

M  Mach  number 

Re  Reynolds  number 

U,  u  Velocity 

u '  Fluctuating  component  of  stream  velocity 

x  Streamwise,  chordwise,  or  longitudinal  coordinate 

y  Coordinate  normal  to  surface 

a  Angle  of  attack 

y  Ratio  of  specific  heats  (1.400  for  air  at  standard  conditions) 

$  Boundary  layer  total  thickness 

$*  Boundary  layer  displacement  thickness 

0  Boundary  layer  momentum  thickness 

X,  A  Sweep  angle 


57 


AEDC-TR-84-23 


v  Kinematic  viscosity 

q  Gas  density 

SUBSCRIPTS 

1,2  Upstream  and  downstream  of  shock,  respectively,  or  before  and  after  tripping, 
respectively 

c  Based  on  chord  length 

e  Local  edge-of-boundary  layer  value 

k  At  top  of  roughness  element  or  at  station  of  trip 

sh  At  shock  location 

t  At  transition 

u,  I  Upper  and  lower  surfaces,  respectively 

x  Based  on  distance  x 

fi  Based  on  5  as  characteristic  length 

5*  Based  on  5* 

6  Based  on  $ 

oo  Free-stream  value 
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